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ABSTRACT 

A theoretical light curve is constructed for the quiescent phase of the recurrent nova U Scorpii in 
order to resolve the existing distance discrepancy between the outbursts (d ~ 6 kpc) and the quiescences 
(d ~ 14 kpc). Our U Sco model consists of a very massive white dwarf (WD), an accretion disk (ACDK) 
with a flaring-up rim, and a lobe- filling, slightly evolved, main-sequence star (MS). The model properly 
includes an accretion luminosity of the WD, a viscous luminosity of the ACDK, a reflection effect of 
the MS and the ACDK irradiated by the WD photosphere. The B light curve is well reproduced 
by a model of 1.37 M Q WD + 1.5 M Q MS (0.8—2.0 M MS is acceptable) with an ACDK having 
a flaring-up rim, and the inclination angle of the orbit i ~ 80°. The calculated color is rather blue 
(B — V ~ 0.0) for a suggested mass accretion rate of 2.5 x 1O~ 7 M yr _1 , thus indicating a large color 
excess of E(B — V) ~ 0.56 with the observational color of B — V = 0.56 in quiescence. Such a large 
color excess corresponds to an absorption of Ay ~ 1.8 and As ~ 2.3, which reduces the distance to 6 — 8 
kpc. This is in good agreement with the distance estimation of 4 — 6 kpc for the latest outburst. Such a 
large intrinsic absorption is very consistent with the recently detected period change of U Sco, which is 
indicating a mass outflow of ~ 3 x lO~'Af yr _1 through the outer Lagrangian points in quiescence. 
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individual (U Scorpii) 
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1. INTRODUCTION 



1989). The residual 40% (1.2 x 10" 6 M Q ) is added to the 



U Scorpii is one of the best observed recurrent no- 
vae, the outbursts of which were recorded in 1863, 1906, 
1936, 1979, 1987, and the latest in 1999. Especially, the 
1999 outburst was well observed from the r ising phase to 
the c o oling phase by many o bservers (e.g. , [Vlunari et al 



1999 



Kahabka et al. 1999; Lcpinc et al. 1999| ) includ 



mg eclipses (Matsumoto, Kato, & Hachisu 2000). Based 
on Matsumoto et al.'s (2000) observation, Hachisu et al. 



(2000) have oonotruotcd a theoretical light curve model 
for the 1999 outburst of U Sco and obtained various prrys 



helium layer of the WD. (5) As a result, the WD can grow 
in mass at an average rate of ~ 1 x 1O~ 7 AT yr" 1 . 

The above physical pictures are exactly the same as pro- 
posed by Hachisu et al. (1999b) as a progenitor system of 
Type la supernovae (SNe la). However, the distance to 
U Sco is still controversial because the direct light-curve 
fitting results in a r elatively short distance of ~ 6 kpc 
(Hachisu et al. 2000), which is incompatible wit h the dis- 
tance of rZ 1 4 kpc at the quiescent phase (e.g . , Webbink 



et al. 1987 



ical parameters of the recurrent nova. Their main results 
are summarized as follows: (1) A direct light-curve fit- 
ting of the 1999 outburst indicates a very massive white 
dwarf (WD) of Af WD = 1.37 ± O.O1M . (2) The en- 
velope mass at the optical maximum is estimated to be 
AM ~ 3 x 1O~ 6 M . (3) Therefore, the mass accretion 
rate of the WD is Af acc ~ 2.5 x 1O~ 7 M yr _1 during the 
quiescent phase between 1987 and 1999. (4) An optically 
thick wind blows from the WD and plays a key role in 
determining t he nova duration becau se it reduces the en- 
velope mass (Kato fc Hachisu 1994). About 60% of the 



mary) . If t 



Warner 1995; Kahabka et al 
ic distance of 



199E , for a sum- 
14 kpc is the case, 



envelope mass is carried away in the wind, wh ich forms an 
expanding shell as observed in T Pyx (e.g., Shara et al 



it could 

be hardly consistent with the results (1) to (5) mentioned 
above. 

Our purpose in this Letter is to construct a light-curve 
model for the quiescent phase and to rectify the distance 
to U Sco. Our numerical method to obtain light curves 
has been described both in Hachisu & Kato (1999) to ex- 
plain the second peak of T CrB outbursts and in Hachisu 
et al. (2000) to reproduce the light curve for the 1999 out- 
burst of U Sco. Therefore, we mention only new parts of 
our numerical method in §2. In §3, by directly fitting our 
theoretical light curve to the observations, we derive the 
distance to U Sco. Discussions follow in §4, especially in 
relation to the recently detected orbital-period change of U 
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Sco and a systemic mass loss through the outer Lagrangian 
points. We also discuss the relation to a progenitor system 
of SNe la. 

2. THEORETICAL LIGHT CURVES 

Our U Sco model is graphically shown in Figure [l| 
Schaefer (1990) and Schaefer & Ringwald (1995) observed 
eclipses of U Sco in the quiescent phase and determined the 
orbital period (P — 1.23056 days) and the cphcmcris (HJD 
2,451,235.777 + 1.23056£) at the epoch of mid-eclipse. 
Thus, the companion is a main-sequence star (MS) which 
expands to fill its Roche lobe after a most part of the cen- 
tral hydrogen is consumed. We call such a star "a slightly 
evolved" MS. The inclination angle of the orbit (i ~ 80°) 
is a parameter for fitting. 

We have assumed that (1) M W d = 1.37M©, (2) the WD 
luminosity of 



Lwd 



1 GMwDM a , 

2 



R 



L 



WD,0, 



WD 



where the first term is the accretion luminosity (e.g., 
Starrfield, Sparks, & Shaviv 1988) and the second term 
Lwu,o is the intrinsic luminosity of the WD, and i?wD = 
O.OO32i? the radius of the 1.37M© WD, and (3) a black- 
body photosphere of the WD. The accretion luminos- 
ity is ~ 1700L© for a suggested mass accretion rate of 
A^acc ~ 2.5 x 10~ 7 Af Q yr^ 1 . Here, we assume Lwd,o = 
because the nuclear luminosity is smaller than the accre- 
tion luminosity for this accretion rate, but we have exam- 
ined other two cases of Lwd.o = 2000 and 4000 L Q and 
found no significant differences in the distance as shown 
below. We do not consider the limb-darkening effect for 
simplicity. 

It is assumed that the companion star is synchronously 
rotating on a circular orbit and its surface fills the in- 
ner critical Roche lobe as shown in Figure [l]. We neglect 
both the limb-darkening effect and the gravity-darkening 
effect of the companion star for simplicity. Here, we as- 
sume a 50% irradiation efficiency of the companion star 
(?7ir,MS = 0.5). We have examined the dependence of the 
distance on the irradiation efficiency (i.e., ?7i ri MS — 0-25 
and 1.0) but found no significant differences in the dis- 
tance as shown below. The non-irradiated photospheric 
temperature Xph.MS of the companion star is a parameter 
for fitting. The mass of the secondary is assumed to be 
M M s = 1.5M Q . 

The size of the accretion disk is a parameter for fitting 
and defined as 

-Rdisk = <xRi, (2) 

where a is a numerical factor indicating the size of the ac- 
cretion disk, and Rl the effective radius of the inner critical 
Roche lobe for the WD component (e.g., Eggleton 1983). 
We also assume that the accretion disk is axisymmetric 
and has a thickness given by 



h = (3R d 



Lsk 



isk 



(3) 



where h is the height of the surface from the equatorial 
plane, w the distance on the equatorial plane from the 
center of the WD, v the power of the surface shape, and 



a numerical factor showing the degree of thickness and 
also a parameter for fitting. We adopt a tu-squared law 
(v — 2) simply to mimic the flaring-up effect of the accre- 
tion disk rim (e.g., Schandl, Meyer-Hofmeister, & Meyer 
1997), and have examined the dependence of the distance 
on the power (y = 1.25 and 3.0) without finding any sig- 
nificant differences as shown below. 

The surface of the accretion disk also absorbs photons 
from the WD photosphere and reemits with a black-body 
spectrum at a local temperature. We assume a 50% irra- 
diatio n efficiency of the co mpanion star, i.e., = 0.5 
(e.g., Schandl et al. 1997). We have examined other two 



cases of ?7i r ,DK = 0.25 and 1.0, and found no significant 
differences in the distance as shown below. The non- 
irradiated temperature of the disk surface is assumed to 
be determined by the viscous heating of the standard ac- 
cretion disk model. Then, the disk surface temperature is 
given by 



<JT, 



ph.disk 



3GM WD M a 
87rn7 3 



where r the distance from the WD center, and cosfl th e 
incident angle of the surface (e.g., schandl et al. 1997). 
The temperature of the disk rim is assumed to be 3000 K. 

3. RESULTS 

Figure || shows the observational points (open circles) by 
Schaefer (1990) together with our calculated B light curve 
(thick solid line) for M acc = 2.5 x 10" 7 Af o yr" 1 . To fit 
our theoretical light curves with the observational points, 
we calculate B light curves by changing the parameters of 
a = 0.5—1.0 by 0.1 step, (3 = 0.05—0.50 by 0.05 step, 
T phi MS = 3500—8000 K by 100 K step, and i = 75—85° 
by 1° step and seek for the best fit model. The best fit 
parameters obtained are shown in Figure ^ (see also Table 

There arc five different contributions to the B-light (Lb) 
in the system: the white dwarf (Lbi), the non-irradiated 
portions of the accretion disk (Lb2) and the donor star 
{Lbs,), and the irradiated portions of the accretion disk 
(Lb a) and the donor star (Lbs)- In order to show each con- 
tribution, we have added two light curves in Figure |2[ that 
is, a non-irradiation case of the ACDK (t^dk = 0, dash- 
dotted), and a non- irradiation case of the MS (?7i r ,MS = 0, 
dashed). The light from the WD is completely blocked by 
the accretion disk rim, thus having no contribution, Lbi = 
0. The depth of the primary eclipse, 1.5 mag, means 
Lb3 = 0.25Lb because the ACDK is completely occulted 
by the MS. The difference of 1 mag between the thick solid 
and dash-dotted lines indicates Lba — O.GOLb- The differ- 
ence of 0.1 mag between the thick solid and dashed lines in- 
dicates Lbs = O.IOLb- Thus, we obtain each contribution: 
Lbi = 0, L B2 = 0.05L B , L B3 = 0.25i B , L Bi = 0.60L B , 
and Lb5 — Q.IOLb- 

Then we calculate the theoretical color index (B — V) c 
for these best fit models. Here, we explain only the case 
of M acc = 2.5 x 10" 7 Af o yr" 1 . By fitting, we obtain the 
apparent distance modulus of tob.o = 16.71, which corre- 
sponds to the distance of d — 22 kpc without absorption 
(Ab = 0). On the other hand, we obtained a rather blue 
color index of (B — V) c = 0.0 outside eclipses. Together 
with the observed color of (B — V) — 0.56 outside eclipses 
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flSchaefer 1990j |Schacfer fe Ringwald 1995| ), we derive a 
color excess of E(B - V) = [B - V) Q - [B - V) c = 0.56 
Here, suffixes c and o represent the theoretically calcu- 
lated and the observational values, respectively. Then, we 
expect an absorption of Ay — 3.1 E(B — V) = 1.8 and 
A B = A v + E(B - V) = 2.3. Thus, we are forced to have 
a rather short distance to U Sco of 7.5 kpc. 

In our case of a — 0.7 and (3 — 0.30, the accretion disk 
is completely occulted at mid-eclipse. The color index of 
(B — V) c = 0.53 at mid-eclipse indicates a spectral type of 



F8 foi the cool component MS, which is in good agreement 
with the spectral type of F8±2 suggested by Johnston & 
Kulkarni (f 992). Hancs (1985) also suggested that a spec- 
tral type nearer F7 is preferred. 

For other mass accretion rates of M acc = (0.1 — 
5.0) x 10 _7 Mq yr _1 , we obtain similar short distances to U 
Sco, as summarized in Table [I], ft should be noted that, al- 
though the luminosity of the model depends on our various 
assumptions of the irradiation efficiencies, the ro-powered 
law of the disk, and the intrinsic luminosity of the WD, 
the derived distance to U Sco itself is almost independent 
of these assumptions, as seen from Table ||. Therefore, the 
relatively short distance to U Sco (~ 6 — 8 kpc) is a rather 
robust conclusion, at least, from the theoretical point of 



4. DISCUSSION 

Matsumoto et al. (2000) observed a few eclipses during 
the 1999 outburst and, for the first time, detected a sig- 



nifica: it period-change of P/P = (-1.7±0.7) x 10~ b yr" 



If we assume the conservative mass transfer, this period 
change requires a mass transfer rate of ;> 10~ 6 M Q yr -1 in 
quiescence. Such a mass transfer for 12 years is too high to 
be compatible with the envelope mass on the white dwarf, 
thus implying a non-conservative mass transfer in U Sco. 

We have estimated the mass transfer rate for a non- 
conservative case by assuming that matter is escaping from 
the outer Lagrangian points and thus the specific angular 
momentum of the escaping matt er is 1.7a 2 f2 or b ( [Sawada 
et al. 1984; Hachisu et al. 1999a), where a is the separa- 



tion and il or b ~~ 2ir/P. Then the mass transfer rate from 
the companion is Mms = {— 5.5 ± 1.5) x 10 _Y M Q yr _i for 
M M s = 0.8 — 2.0 M Q under the assumption that the WD 
receives matter at a rate of M acc = 2.5 x 10~ 7 M Q yr -1 . 
The residual (~ 3 x 10~ 7 M Q yr -1 ), which is escaping from 
the system, forms an excretion disk outside the orbit of the 
binary. Such an extended excretion disk/torus may cause 
a large color excess of E(B — V) = 0.56. 

Kahabka et al. (1999) reported the hydrogen column 
density of (3.1 — 4.8) x 1 21 mi' 2 , which is much larger 



than the Galactic absorption in the direction of U Sco al. 1984 ) 



(1.4xl0 21 cm -2 , Dickey & Lockman 1990), indicating a 
substantial intrinsic absorption. It should also be noted 
here that Barlow et al. (1981) estimated the absorption 
toward U Sco by three ways: (1) the Galactic absorption in 
the direction of U Sco, E(B-V) ~ 0.24 and A v ~ 0.7, (2) 
the line ratio of He II during the 1979 outburst (t ~ 12 days 
after maximum), E(B — V) ~ 0.2 and Ay ~ 0.6, and (3) 
the Balmer line ratio during the 1979 outburst (t ~ 33 — 
34 days after maximum), E(B — V)~ 0.35 and Ay ~ 1.1. 
The last one is significantly larger than the other two es- 



timates. They suggested the breakdown of their case B 
approximation in high density regions. However, we may 
point out another possibility that the systemic mass out- 
flow from the binary system has already begun at t ~ 33 
days and, as a result, an intrinsic absorption is gradually 
increasing. 

The mass of the companion star can be constrained 
from the mass transfer rate. Such a high transfer rate 
as Mms ~ 5.5 x 10 _7 A/ Q yr -1 str ongly indicates a ther- 



mally unstable mass transfer (e.g., van den Heuvel et al 



1992), which is realized when the mass ratio is larger than 
1.0 — 1.1, i.e., q — Mms/Mwd > 1 -0 — 1-1 for zero-age 
main-sequence stars (Webbink 1985). This may pose a 
requirement Mms > 1.4M©. We estimate the most likely 
companion mass of 1.4 — 1.6 M from equation (11) in 
Hachisu et al. (1999b). 

If the distance to U Sco is ~ 6.0 — 8.0 kpc, it is located ~ 
2.3 — 3.0 kpc above the Galactic plane (b = 22°). The zero- 
age masses of the progenitor system to U Sco are rather 
massive (e.g., 8.0 M Q + 2.5 M from Hachisu et al. 1999b) 
and it is unlikely that such massive stars were born in the 
halo. Some normal B-type main-sequence stars have been 
found in the halo (e.g., PG0009+036 is loca ted ~ 5 kpc 
below the Galactic disk, Schmidt et al. 1996), which were 
ejected from the Galactic disk because of their relatively 
high moving velocities ^100 — 200 km s" 1 . The radial ve- 
locity of U Sco is not known but it is suggested that the 
7- velocit y is ^50 — 100 km s -1 from t he absorption line ve- 
lociti es (Johnston & Kulkarni 1992; Bchaefer fc Ringwaldj 



1995). If so, it seems likely that U Sco was ejected from 
the Galactic disk with a vertical velocity faster than ~ 20 
km s _1 and has reached at the present place within the 
main-sequence lifetimes of a ~ 3.0M Q star (~ 3.5 x 10 8 

yr). 

Now, we can understand the current evolutionary sta- 
tus and a further evolution of U Sco system. The white 
dwarf has a mass 1.37 ± 0.01M Q . It is very likely that the 
WD has reached such a large mass by mass accretion. In 
fact the WD is currently increasing the mass of the helium 
layer at a rate of Mh ~ 1-0 X lO~ 7 A/ yr -1 ( |Hachisu~e"t 
al. 2000). We then predict that the WD will evolve as 
follows. When the mass of the helium layer reaches a crit- 
ical mass after many cycles of recurrent nova outbursts, 
a helium shell flash will occur. Its strength is as weak as 
those of AGB s tars because of the high mass accretion rate 
(Nomoto 1982). A part of the helium layer will be blown 
off in the wind, but virtually all of the helium layer will 
be bur nt into carbon-oxygen a nd accumulates in the white 
dwarf (Kato & Hachisu 1999). Therefore, the WD mass 
can grow until an SN la explosion is triggered (Nomoto et 
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Table 1 
u sco quiescent phase? 



-^acc 

(Mq yr- 1 ) 


a 




Tph.MS 

(K) 


"is,0 


outside eclipse 
(B - V)c 


midcclipse 
(B - V) c 


outside eclipse 
E(B - V) 


Ay 


A B 


d 

(kpc) 


5.0X10" 7 


0.7 


0.30 


5900 


17.18 


-0.08 


0.45 (F5) 


0.64 


2.01 


2.65 


8.0 


2.5xl0- 7 


0.7 


0.30 


5500 


16.71 


+0.00 


0.53 (F8) 


0.56 


1.76 


2.32 


7.5 


l.Oxltr 7 


0.7 


0.25 


5000 


16.02 


+0.12 


0.66 (G4) 


0.44 


1.38 


1.82 


6.9 


5.0X10" 8 


0.7 


0.25 


4600 


15.45 


+0.24 


0.78 (G9) 


0.32 


1.01 


1.33 


6.7 


2.5xl0- 8 


0.7 


0.25 


4200 


14.72 


+0.37 


0.91 (K2) 


0.19 


0.60 


0.79 


6.1 


LOxlO" 8 


0.7 


0.25 


3700 


13.58 


+0.58 


1.13 (K5) 








5.2 



inclination angle i = 80° for all cases 



Table 2 

Model dependence of the distance. 



modcl a 


Tph.MS 

(K) 


i rn B ,o 


outside eclipse 
(B - V) c 


midcclipse 

(B - V) c 


outside cclipst 
E(B - V) 


3 A V 


A B 


d 

(kpc) 


L WD ,o = 2000L Q 


6100 


17.31 


-0.09 


0.41 (F3) 


0.65 


2.04 


2.60 


8.4 


Iwd,o = 4000L Q 


6400 


17.59 


-0.13 


0.36 (Fl) 


0.69 


2.17 


2.87 


8.8 


»?ir,MS = 10 


5500 


16.76 


+0.00 


0.53 (F8) 


0.56 


1.76 


2.32 


7.7 


riii,MS = °- 25 


5500 


16.66 


+0.00 


0.53 (F8) 


0.56 


1.76 


2.32 


7.4 


Vii,DK = 10 


6000 


17.16 


-0.06 


0.43 (F4) 


0.62 


1.95 


2.57 


8.3 


Vii,DK = 0.25 


5300 


16.39 


+0.08 


0.58 (F9) 


0.48 


1.51 


1.99 


7.6 


v = 3.0 


5600 


16.90 


-0.04 


0.51 (F7) 


0.60 


1.88 


2.48 


7.6 


v = 1.25 


5300 


16.39 


+0.09 


0.58 (F9) 


0.47 


1.47 


1.95 


7.7 


H = 80°, Mace = 2.5 x 10- 
assumcd, otherwise specified. 


7 M Q yr- 1 , 


a = 0.7, = 0.30, 


V = 2, »?ir,MS 


= 0.5, »?ir,DK 


= 0.5, and Lwd.o 


= arc 
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Fig. 1. — Configuration of our U Sco model at quiescent phase. The cool component (right figure) is a slightly evolved MS (1.5Mq) filling 
up its inner critical Roche lobe. Only the north and south polar areas of the secondary are slightly heated up by the hot component (1.37 Mq 
WD, left figure) because a large part of the light from the hot component is blocked by the flaring-up edge of the accretion disk. Both the 
hot component and the central part of the accretion disk are not seen from the Earth because they are blocked by the flaring-up rim. Here 
the separation is a = 6.87Rq, the effective radii of the inner critical Roche lobes are R^ = 2.55i?0, and R^ = R2 = 2.66R0, for the primary 
WD and the secondary MS, respectively. 
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Fig. 2. — Top: theoretical B light curves (thick solid line) and Bottom: theoretical B — V light curves (thin solid line) plotted against 
the binary phase (two phases from —1.0 to 1.0) together with the observational points (open circles in the B light curve represent data taken 
from Schaefer 1990, but no data points in the B — V light curve). The model is a binary system of 1.37M0 WD + 1.5Mq MS. The other 
parameters are printed in the figure. Two light curves are added in order to clarify each contribution of the light, for the case of no irradiation 
of the ACDK, t)i r] DK = (dash-dotted), and for the case of no irradiation of the MS, r/j ri MS = (dashed). 



